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On the basis of first-principles computational approaches, we present a new method to drive 
zigzag graphene nanoribons (ZGNRs) into the half-metallic state using a ferroelectric material, 
poly(vinylidene fluoride) (PVDF). Owing to strong dipole moments of PVDFs, the ground state of 
the ZGNR becomes half-metallic when a critical coverage of PVDFs is achieved on the ZGNR. Since 
ferroelectric polymers are physisorbed, the direction of the dipole field in PVDFs can be rotated 
by relatively small external electric fields, and the switching between half-metallic and insulating 
states may be achieved. Our results suggest that, without excessively large external gate electric 
fields, half-metallic states of ZGNRs are realizable through the deposition of ferroelectric polymers 
and their electronic and magnetic properties are controllable via noninvasive mutual interactions. 



Since graphene, a single layer of graphite, was shown to 
be exfoliated on the sihcon oxide surface^, there has been 
many studies on its various interesting physical, chemi- 
cal, and mechanical properties. Specially, its unique elec- 
tronic characteristics have spurred researchers to envis- 
age nanoelectronic circuits composed of carbon atoms 
only^i^. To achieve nanoelectronic devices based on 
graphene, it should be cut into small pieces which in- 
evitably have edges at the boundaries. Together with the 
finite size effects^i^, nanoscale graphene fragments with 
these edges exhibit many different properties compared 
with those of two-dimensional graphene. Thus, it is im- 
portant to understand the physics of the one-dimensional 
graphene nanostructure, that is, graphene nanoribbon 
(GNR)^i^. Among various properties of the ribbon stud- 
ied^""—, the magnetism of the zigzag graphene nanorib- 
bon (ZGNR) is especially notableiir— . The magnetism 
arises at the zigzag-shaped edges because of the localized 
TT-orbitals of carbon atoms at the edges^iii"— . Consid- 
ering a quite long spin coherence length in graphene^^, 
the magnetism at the edges may be useful in future spin- 
tronics applications^. 

One of the most interesting properties regarding the 
application of the magnetism in ZGNRsi^"— is the 
electric-field-induced half-metallicity^^. Under the trans- 
verse electric field, ZGNRs show the half-metallicityii 
which originates from a unique interplay between the in- 
teredge antiferromagnetic ordering and the relative po- 
tential shift between two edges. Despite its novel fea- 
ture, the required transverse electric field for the half- 
metallicity is too strong to be obtained easily in ex- 
perimentsii. Many researchers have suggested diverse 
methods, such as introducing quantum dots with the 
zigzag edgesi^iii, functionalizing edges of ZGNRs^®"^°, 
and boron-nitrogen (B-N) substituting for designated 
carbon atoms^^'^^. Some of these methods are success- 
ful in showing half-metallicity without external electric 
fields, though impractical because fine control of posi- 
tions of functional groups and B-N substitution is in- 



evitable. Moreover, the half-metallic states are hardly 
controllable with those proposed methods. 

In this paper, we propose a method to drive the ZGNR 
into a half-metallic state by depositing a ferroelectric ma- 
terial, poly(vinyhdene fiuoride) (PVDF)^^ on the ZGNR. 
The PVDF, one of the well-known ferroelectric materials, 
shows strong ferroelectricity and high crystallinity^''"^''. 

The atomic model for the PVDF is in Figure la. In 
general, two surfaces of the ferroelectric thin film, cut 
across the dipole moment aligned to a certain direc- 
tion, cause a potential difference due to the dipole field. 
Likewise, several aligned PVDFs form the electric field 
which is induced by the dipole moments, and even a 
few layers show a substantial ferroelectricity. Here, we 
show that PVDFs can weakly bind to the surface of 
the ZGNR, with dipole direction parallel to the surface, 
thereby achieving the half-metallicity without any elec- 
tric field. Because of their weak binding nature, the di- 
rection of induced electric fields is controllable, and the 
switching between half-metallic states and antiferromag- 
netic insulating states are made possible by external gate 
fields. We also note other experiment studies exploiting 
the interplay beetween ferroelectricity of the PVDF and 
graphene^S^. 



Results and Discussion 

We have considered the n-ZGNR that has n units of 
the C-C pair in the x-direction per unit cell and repeats 
along the y-axis with the periodicity of 2.46 A. The unit 
cell of the PVDF also repeats along the y-axis, but with 
a different equilibrium unit cell length of 2.54 A which 
is 3.15% longer than that of the ZGNR (Figure la). For 
convenience, the unit cell size of the PVDF is set to be 
reduced in the y-direction so that it matches the length 
of the ZGNR (2.46 A) in studying the joint system (Fig- 
ure lb,c). The unit cell in this system is indicated in 
Figure Ic. We confirm that our atomic models with the 
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FIG. 1: (a) Ball-and-stick model for Poly (vinylidene fluo- 
ride) (PVDF). The arrows indicate the direction of the dipole 
moment (from the fluorine side to the hydrogen side), (b) 
Stable configuration of a PVDF chain deposited on the 8- 
ZGNR whose edges are passivated by hydrogen atoms. The 
PVDF chain and the ZGNR are set parallel along the y- 
axis. The PVDF lies flat on the ZGNR, the height from 
the ZGNR to the lowest atom of the PVDF is 2.5 A, and 
the PVDF sits on the center of the ZGNR. (c) Top view of 
the same configuration as in (b). (d) The isosurface of the 
charge density difference between the up-spin and the down- 
spin [p^{r) — Pi{r)) when the PVDF is deposited on the 8- 
ZGNR. The blue (red) region denotes the + sign (— sign) 
with the isovalue of 6.7 x 10"^|etA~^. 



slightly compressed PVDF on the ZGNR give essentially 
the same results as the ones obtained from a fully relaxed 
geometry (without compression) in a very large commen- 
surate unit-cell for the PVDF and the ZGNR. We will 
discuss this issue later. The vacuum along the x-axis is 
set to be longer than 70 A in order to avoid the spurious 
dipole interactions between the adjacent PVDFs in re- 
peated supercells. All properties remain the same when 
the supercell size is doubled, and the dipole correction 
is negligible for our purpose. In determining the most 
stable position of the PVDF on the ZGNR, three param- 
eters have to be optimized simultaneously: directions of 
the dipole moment of the PVDF, heights of the PVDF 
from the ZGNR plane, and distances between the PVDF 
and the edge of the ZGNR along the x-axis. From a com- 
binatorial search, the structure is found to be stabilized 
when the dipole orientation is parallel to the surface of 
the ZGNR, the height of the lowest atom of the PVDF 
from the ZGNR is 2.5 A, and the PVDF lies on the mid- 
dle position of the ZGNR (Figure lb,c). 

From theoretical*' and experimenta P'^i^^" — studies, 
the dipole orientation of the PVDF is known to depend 
highly on the kind of substrate. When the substrate is a 
conductor, the image dipole is induced by the dipole of 
the PVDF on conducting substrates and the interaction 
energy between two dipoles should be taken into account. 
If we regard the PVDF as a single dipole, the energy of 
perpendicular dipoles is ~2p'^/r^ from the dipole-dipole 
interaction energy equation'^*', where p and r denote the 
magnitude of the dipole moment and the distance be- 



tween the dipole of the PVDF and its image dipole, re- 
spectively. On the other hand, the energy of parallel 
dipoles is —p-^/r^. Thus, on the conducting substrate, 
the perpendicular dipole direction of the PVDF is en- 
ergetically favorcd'^"-'^^. On the contrary, in the case of 
an insulating or semiconducting substrate, for example, 
KBr— or KCSS,, the dipole is aligned parallel to the sur- 
face plane because the image dipole is negligible^° . From 
the calculations, we find that the PVDF favors lying flat 
to the ZGNR plane since the ZGNR is semiconducting 
with its low screening capability. 

The distance from the ZGNR to the lowest atom of 
the PVDF is found to be 2.5 A, and it is 3.5 A from the 
ZGNR to the C-C bond of the PVDF (Figure lb). From 
the fully relaxed geometries, we can infer that the PVDF 
weakly binds to the ZGNR and the PVDF affects the 
electronic structure of the ZGNR not through chemical 
bonding between the PVDF and the ZGNR, but through 
the electronic potential induced by the strong dipole mo- 
ment. The structure shown in Figure Ic, where the C-C 
bonds of the PVDF and those of the ZGNR cross each 
other at the middle of the ZGNR, has been found to 
a ground state. Even if the PVDF deviates its ground 
state geometry slightly, we find that no essential varia- 
tion occurs in the electronic structures due to the weak 
bonding between the PVDF and the ZGNR. It is also no- 
ticeable that, with the PVDF, the antiferromagnetic spin 
configuration of the ZGNR is still more stable than the 
ferromagnetic one. For a single PVDF on the 8-ZGNR, 
the interedge antiferromagnetic ordering is more stable 
by 2.1 meV per edge atom. The isosurface of the charge 
density difference between the up-spin and the down-spin 
(/9-|-(r) — Pi{r)) in this case is presented in Figure Id. 

Now, we analyze the band structure of this joint sys- 
tem, a single PVDF on the 8-ZGNR. The pristine 8- 
ZGNR shows a small band gap of 0.3 eV for both spins. 
With one PVDF deposited on the 8-ZGNR, the band 
structure of the ZGNR is changed as if it were under ex- 
ternal transverse electric fields (Figure 2). The electric 
potential energy is raised at the fluorine side and lowered 
at the hydrogen side of the PVDF. Since the dipole field 
decays slowly in space, the potential drop induced by this 
dipole field takes place over a large area which can cover 
both edges of the ZGNR. The energy eigenvalues of the 
ZGNR are also modified by this potential variation so 
that the band gap of the, say, up-spin state is reduced 
to 0.18 eV, and that of the down-spin state is increased 
up to 0.41 eV (Figure 2). We can estimate the potential 
profile generated by a single PVDF chain on the ZGNR 
by averaging the potential energies in the y-axis with the 
position of the ZGNR fixed in z direction. Its drop across 
two edges of the PVDF is 1.2 V and that of the 8-ZGNR 
is 0.8 V. According to previous works^^, such a band gap 
change corresponds to the case of homogeneous electric 
field of 0.05 V/A applied to the 8-ZGNR, that is, the to- 
tal potential difference between two edges of the 8-ZGNR 
is 0.89 V ( = 0.05 V/A x 17.85 A, where 17.85 A is the 
length of the 8-ZGNR), which agrees reasonably well with 
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FIG. 2: The electronic band structures near the Fermi level. 
From left to right, the number of the PVDF chains is increased 
from one to four. Due to the potential drop by PVDFs, the 
band gap of the up-spin decreases, and that of the down-spin 
increases. The red (blue) line denotes the band structure of 
the up (down) -spin as a function of k. (They overlap in most 
places.) The inset displays the band structure in the range of 
\E - Ef\ < 100 meV and O.Ttt < ka < ty. In aU figures, the 
Fermi level (Ep) is set to zero. 
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FIG. 3: Dependence of half-metallicity on different sizes 
of the ZGNR and the PVDF. A^(red) denotes the direct 
bandgap of the up-spin, and A^(blue) the gap of the down- 
spin as a function of the coverage of PVDFs for the 8- 
ZGNR (filled circles), 12-ZGNR (open circles), and 16-ZGNR 
(squres). When the coverage reaches 16.0/to x 100(%), where 
w is the width of the ZGNR in angstrom, the ZGNR becomes 
half- metallic. 



the present calculation of 0.8 eV. Thus, the change of the 
bandgap mainly originates from the potential difference 
between two edges of the ZGNR. When the electric field 
is applied, the energy eigenvalues of the states localized 
on the left side of ZGNRs are pulled up, and those on 
the right side are pushed down^i. Then, the energy gap 
of the up-spin state decreases, and that of the down-spin 
increases. We conclude that the energy eigenvalues of 
ZGNRs are modified because of interplay between the 
antiferromagnetic interedge ordering and the dipole field 
of the PVDF. 

Before proceeding further, we check the lattice mis- 
match problem between the PVDF and the ZGNR in 
the unit cell. As we mentioned before, the lattice of 
the PVDF is reduced to that of the ZGNR in the y- 
axis. To verify the effect of the compressed PVDF on 
the ZGNR, same calculations are performed for a large 
supercell amounting to 32 unit-cells of the ZGNR and 31 
unit-cells of the PVDF, in which the bonds experience al- 
most no strain. Very little change of the band structure 
is observed with a different width of the PVDF along the 
periodic direction. Hence, the role of the PVDF is only 
to generate the electric field by its dipole moment, and 
the position of the PVDF or the reduced cell size does not 
appreciably infiuence the electronic and magnetic prop- 
erties of the system. In addition, a typical domain size of 
the ferroelectric thin film of the PVDF is around 200 nm 
so that, if considering a nanoscale channel device geom- 
etry, the effect of phase changes along the PVDF would 
not be serious in the present casci^^ 

Next, we increase the number of PVDFs deposited on 
the 8-ZGNR one by one (Figure 2). Dipole moments of 
PVDFs are all aligned in the same direction. The total 
dipole moment of PVDFs is increased, and the poten- 
tial drop across the ZGNR becomes 1.26, 2.0, and 3.0 
V if the number of PVDFs is 2, 3, and 4, respectively, 
in the 8-ZGNR case. (The local potential drop between 



two edges of PVDFs increases to 2.1, 3.4, and 5.0 V as 
increasing the number of PVDFs.) The band gap of the 
up-spin state is decreased while that of the down-spin 
state increases as the number of PVDFs increase. Even- 
tually, if four PVDFs are coated on the 8-ZGNR, the 
system turns to a half-metal. The up-spin states cross 
the Fermi level while the down-spin states still have a 
band gap of 0.37 eV. Four PVDFs generate enough of 
a potential difference (3.0 eV) between two edges of the 
ZGNR to give rise to the half-metallicity of the ZGNR. 
With low defect concentrations at the edges, the half- 
metallic nature of ZGNRs is already shown to be robust 
under the applied external transverse electric field in the 
previous work»ii Since PVDFs here play the same role 
as do external fields in ref 11, we note that the induced 
half-metallicity of the ZGNRs will occur as long as weak 
edge defects or imperfections do not alter the magnetic 
ordering at the edges. 

Our calculations indicate that four PVDFs create a po- 
tential drop of 3.0 V across the n-ZGNR irrespective of 
the width (n value) of ZGNRs (for n > 8, tested up to 
n = 32). The n-ZGNR becomes half-metallic when the 
potential drop across the system generated by PVDFs 
reaches 3.0 V, which has been suggested to be the criti- 
cal value for half-metallicity in the previous work^^. In 
terms of the percentage coverage (i.e., the area of covered 
PVDFs divided by the area of the ZGNR), the critical 
value for the transition to half-metal can be calculated by 
the width of four PVDFs (16.0 A) divided by the width 
of the ZGNR (w), 16.0/u; x 100(%) {w is in angstrom). 
This can be regarded as a scaling behavior in the PVDF 
coverage and the width of ZGNRs as follows, similar to 
the previous study on scaling in the strength of transverse 
electric fields and the width of the ZGNR" (Figure 3). 

Finally, another advantage of our proposal compared 
with others lies in the fact that the half-metallic state can 
return to the insulating state by changing of the dipole 
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orientation of PVDFs. If the dipole direction of PVDFs is 
set perpendicular to the ZGNR plane by a perpendicular 
external electric field as in the usual experiments with 
back gates, then the potential drop across the ZGNR 
vanishes and the half-metallicity disappears as well. We 
have calculated variations of the total energy by rotating 
a single PVDF on ZGNRs. We have found that, without 
a perpendicular electric field, the perpendicular dipole 
direction of the PVDF to the ZGNR is also at the local 
energy minimum (or quasi-stable configuration), 30 meV 
higher than the ground state (the parallel dipole direc- 
tion). Upon application of perpendicular electric fields, 
the energy differences between parallel and perpendicular 
dipole configurations are decreased and eventually, with 
an electric field of '^0.07 V/A, the perpendicular dipole 
direction becomes the ground state. An energy barrier 
for 90° rotation of the PVDF on the ZGNR, that is, par- 
allel to perpendicular direction change, is 96 meV per 
unit cell without an electric field and decreases to ^11 
meV with an electric field of 0.07 V/A. Therefore, it is 
practically possible to switch between half-metallic and 
insulating states by applying a perpendicular gate field 
greater than 0.07 V/A. 



Conclusion 

We have proposed a method to construct the half- 
metallic ZGNRs by depositing the ferroelectric material 
PVDFs on them. The PVDFs generate the electrostatic 
potential on the ZGNR because of their strong dipole mo- 
ments. When the number of PVDFs coated on ZGNRs 
increases beyond a critical value, the system becomes be 
half-metallic. The dipole direction of PVDFs is change- 
able by an external electric field, and switching between 
half-metallic and insulating states in ZGNRs is made pos- 
sible. 



Theoretical Methods 

We performed the first-principle calculations based 
on the density-functional theory (DFT) within the lo- 
cal spin density approximation (LSDA)^ and spin po- 
larized general gradient approximation (GGA) with 
the Perdew-Burke-Ernzcrhof (PBE) functional2^, respec- 
tively, by using the SIESTA package^. The stan- 
dard norm-conserving TrouUier-Martins pseudopoten- 
tials'^^ were employed, and split-valence double-C plus 
polarization basis'^^ was used. We have chosen a 400 Ry 
energy cutoff for a real space mesh size and 96 /c-points, 
uniformly distributed in the ID Brillouin zone. All edges 
were saturated with hydrogen atoms and relaxed by con- 
jugate gradient minimization until the maximum force 
was less than 0.04 eV/A. To overcome an intrinsic error 
due to the pseudoatomic orbital basis set in the weakly 
binding system, the basis set superposition error (ESSE) 
was corrected using a counterpoise procedure3S.i22.. In 



calculations with the PBE functional, we added a pair- 
wise interatomic CqR^^ term (-Evdw) to the PBE-DFT 
energy in order to include the van der Waals (vdW) in- 
teraction^"—. To obtain accurate vdW energies, the ef- 
fective Cg coefficients of each atom in the systems were 
obtained by using a recent theoretical method exploit- 
ing ground state electron density from A. Tkatchencko 
and M. Scheffler"*^. From the calculation including PBE 
and vdW corrections, we found that there is no signifi- 
cant difference between atomic and electronic structures 
based on LDA and those based on PBE and vdW cor- 
rections. The comprehensive methods and comparisons 
are in the Supporting Information. Hence, in the article, 
all binding energies, relative distances, and orientations 
between PVDFs and ZGNRs and band structures in the 
paper are obtained within the LDA calculations. 
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SUPPLEMENTARY INFORMATION 

It is well known that the long range interaction tail^ 
is not included in local-density approximation (LDA)^ 
and generalized gradient approximation (GGA)^i^ to 
exchange-correlation energy functional in density func- 
tional theory (DFT). Because either LDA or GGA are 
based on the uniform electron gas system, they fail to 
describe weak interactions between electrons in sparse 
materials such as layered systems, liquid crystals, poly- 
mers, proteins, and biomolecular surfaces^. Neverthe- 
less, LDA and GGA calculations provide upper and lower 
limits of bond lengths, binding energies, and energy band 
gaps. They can be useful for reference data sets. More- 
over, many excellent methods^"— i^"— which can de- 
scribe the long range interaction or van der Waals (vdW) 
forces within DFT calculations are proposed. 

In our main article. Poly (Vinyliden Fluoride) (PVDF) 
can be deposited on the zigzag graphene nanorib- 
bon (ZGNR) with the parallel dipole direction to the 
graphene plane. From binding energies and bond lengths 
between them from LDA and GGA approximations, 
we infer that PVDFs weakly bind to ZGNR, therefore 
the vdW correction should be included. In this sup- 
porting information, we present binding energies and 
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FIG. 4: Binding energies as increasing the distance between 
8-ZGNR and PVDF from 2.0 A to 4.6 A. We carried out 
LDA (dotted line), GGA (dashed hne) and vdW corrected 
GGA, it is called GGA-D (solid line), (a) With parallel dipole 
direction of PVDF to 8-ZGNR and (b) with perpendicular 
dipole direction of PVDf to 8-ZGNR. 



bond lengths from LDA and GGA calculations when one 
PVDF is deposited on the 8-ZGNR. The vdW correction 
is considered by adding a pairwise interatomic CqR^^ 
term (Eydw) to the DFT energjj^"— i^. Finally, com- 
parisons with LDA and GGA for the ground state under 
the external electric fields are also discussed. 

We performed LDA and GGA calculations when de- 
positing one PVDF on 8-ZGNR as changing the dis- 
tance between the two from 2.0 A to 4.6 A and the 
dipole direction of PVDF, parallel and perpendicular to 
8-ZGNR. In order to reduce an intrinsic error due to the 
small number of basis sets, basis set superposition error 
(BSSE), we consider corrections for the error by intro- 
ducing ghost orbital basis sets^i^. We compare binding 
energies at each distance with LDA and GGA, respec- 
tively. Binding energies are energy difference between to- 
tal system (PVDF deposited on 8-ZGNR) and separated 
system (8-ZGNR, PVDF). They are calculated by fol- 
lowing formula, Eund.cp = £^(8-ZGNR-f PVDF) - £;(8- 
ZGNR, ghost atoms of PVDF) - E{VVT)¥, ghost atoms 
of 8-ZGNR), CP means 'counterpoise' correction^, and 
ghost atoms have only pseudo atomic orbitals without 
potential'^^. 

In Fig. 1(a), dotted (dashed) line represents binding 



Of 



H1 



F6C 1„H2 
018 C16 014 



■jp 

C19 c'l? C15 C13 




3 3.5 
distance (A) 

-^free 



and 8-ZGNR as increasing distance between them, (a) The 
ratio for all atoms in PVDF, and (b) for some atoms in 8- 
ZGNR which are close to PVDF. The number index of atoms 
in each figure is denoted in the atomic model above the panels. 



energies calculated by LDA (GGA) with parallel dipole 
of PVDF to 8-ZGNR. The binding energy of PVDF on 
8-ZGNR from LDA (GGA) is 145 meV (30 meV) and the 
distance between them is 2.5 A (3.2 A) (parallel dipole 
direction to the surface, see inset of Fig. 1(a)). These 
values agree well with tendency of LDA and GGA cal- 
culations, i.e., overbinding in LDA and under binding in 
GGA. When dipole direction of PVDF is perpendicular 
to 8-ZGNR, binding energies are also plotted in Fig. 1(b) 
under LDA and GGA calculations. The binding energy 
is 77 meV and 8 meV, and the distance between PVDF 
and 8-ZGNR is 3.0 A and 3.6 A for LDA and GGA cal- 
culations, respectively. From those results, we confirm 
that parallel dipole of PVDF is more stable than perpen- 
dicular dipole of PVDF on 8-ZGNR irrespective of the 
approximations. 

In order to consider the vdW interaction between 
PVDF and 8-ZGNR, we added vdW correction energy 
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TABLE I: Free-atom Ce coefficients (hartree • bohr^) from 
Chu*^, vdW radius Rfree (A) from Bondi^^ and static polar- 
izability from Chu^ and Yai*^. 

Evdw to the DFT total energy {Edft-d), 

Edft~d — Eks-dft + E^dw, (1) 

where Eks-dft is the usual self-consistent Kohn-Sham 
energy. The vdW energy (-E^dvy) is given by 

EvdW = —7^ ^ fdmp{RAB, Raj R%:)C&AbRaB- (2) 
^ A,B 

Here, Rab is the distance between atoms A and B. Cqab 
is the Cg coefficient between atom A and B which can 
be written as Cqab = ^^s^-^g^g^ where CeAA(BB) 

—^-CeAA^ §-C6BB 

"A "A 

and ce^A(B) ^^'^ homonuclear Cg coefficient and static po- 
larizability of atom A(B). R\ and R^ are the vdW radii. 
The Rab singularity at small distances is eliminated by 
the short-ranged damping function fdmp{RAB , Ra^ R'b)- 

fdmp {Rab, Rab ) = — r r,„ 7TT ' (3) 

1 + exp -d - 1 

where Rab ~ R^a + -^s- ^'^'^ ^'"^ ^^^^ parameters. 
To calculate E^^w, we use a method proposed by A. 
Tkatchenko and M. Scheffler^^. In the method, we can 
calculate the effective parameters, C\^^ and E^ij^ corre- 
sponding to G^AB and i?° respectively from the ground 
state electronic density. This improves the description of 
weakly bonded systems compared with empirical ones^. 

We added this vdW correction to GGA calculations 
with Perdew-Burke-Ernzerhof (PBE)^^ functionals using 
well-known free parameters, d = 20 and sfi — 0.94^, 
Cl^'"^ , R'free' ^^d a^, shown in Table L The summation 
in Eq. (2) is done for atoms within 500 A. 

The binding energy including the ESSE correc- 
tion28, Eund-D = E{8-ZG1^R+PVDF) DFT-D - E{8- 
ZGNR) DFT- D " E(PYDF)dft-d, is plotted as solid 



line in Fig. 1(a) and (b). It gives the binding energy 
of 129 meV and the distance of 2.8 A in (a) and 73 meV 
3.2 A in (b). They are between LDA and GGA limits and 
quite close to the LDA results. From all results, never- 
theless vdW correction changes binding energy and bond 
length between PVDF and 8-ZGNR, it is not a crucial 
part in determining the bonding nature and electronic 
structures. 

Additionally, the effective Cg coefficients, Cq'^^ in the 
system decrease from those of free atoms and changes 
when system configurations are altered. If we change 
the height between PVDF and ZGNR, the C^^^ between 
the nearest pair atoms belong to each material decreases 
significantly while those for well separated pair atoms de- 
creases slightly [Fig. 2] . Because different hybridizations 
states for atoms in each system give various Cg coefficien- 
sts, all Cg values are changed, respectively. The ratio of 
Cg^^ to Cg^''"" is plotted in the Fig. 2(a) and (b). The 
ratio for all PVDF atoms are shown in Fig. 2(a) and the 
ratio for some 8-ZGNR atoms which are close to PVDF 
are presented in Fig. 2(b). The ratios are in the range of 
0.37 to 0.6 for Hydrogen atoms, 0.64 to 0.75 for Carbon 
atoms, and about 0.81 for Fluorine atoms. Especially, 
the closest atom between PVDF and 8-ZGNR, HI (Fig. 
2(a)), has large variation of the Cg ratio, and carbon 
atoms at both edges of 8-ZGNR have constant ratio of 
0.64. 

Finally, the vdW correction is also performed under the 
external electric fields. We found that the LDA results 
are still trustful. Based on electrostatic forces, the dipole 
of PVDF tends to align with parallel direction to external 
electric fields. With the electric fields of 0.1 V/A whose 
direction is perpendicular to ZGNR surface, the perpen- 
dicular PVDF to ZGNR surface is more stable than the 
parallel one. The energy differece between perpendicular 
and parrallel direction of PVDF is 14 meV and 8 meV in 
LDA and vdW corrected GGA, respectively. 

In conclusion, we carried out LDA, GGA, and vdW 
corrected GGA calculations for a single PVDF on 8- 
ZGNR. Our vdW corrected GGA-PBE calculatons show 
that binding energies, distances between PVDF and 
ZGNR, and the response under the external electric fields 
lead to the same conclusion to ones from LDA calcula- 
tions with very small quantitative differences. 
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